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Abstract 

We study the finite size scaling behaviour of the specific heat of thin films 
in the neighbourhood of the A-transition. To this end we have simulated the 
improved two-component (/>^ model on the simple cubic lattice. We employ 
free boundary conditions in the short direction to mimic the vanishing order 
parameter at the boundaries of a '^He film. Most of our simulations are 
performed for the thicknesses Lq = 8, 16 and 32 of the film. It turns out that 
one has to take into account corrections oc Lq to obtain a good collapse of 
the finite size scaling functions obtained from different Lq. Our results are 
compared with those obtained from experiments on thin films of ^He near the 
A-transition, from field theory and from previous Monte Carlo simulations. 



1 Introduction 



In the neighbourhood of a second order phase transition the behaviour of various 
quantities is governed by power laws. For example the correlation length diverges 
as 



where t = (T—Tc) /T^ is the reduced temperature, ^o,+ and ^o,- are the amplitudes in 
the high and the low temperature phase, respectively, and z/ is the critical exponent 
of the correlation length. Tc is the critical temperature, where the phase transition 
occurs. The specific heat behaves as 



where and are the amplitudes in the high and the low temperature phase, re- 
spectively, and B is an analytic background, which has to be taken into account here, 
since the critical exponent a of the specific heat is negative for the three-dimensional 
XY universality class. A universality class is characterised by the dimension of the 
system, the range of the interaction and the symmetry of the order parameter. Crit- 
ical phenomena can be understood in the framework of the Renormalization Group 
(RG). For reviews on critical phenomena and the Renormalization Group see e.g. 
[B [21 H] . The XY universality class in three dimensions with short range interac- 
tions is of particular interest, since the A-transition of ^He shares this universality 
class. Results for critical exponents and amplitude ratios such as Aj^/A- obtained 
at the A-transition, are more precise than those from experiments on other systems. 

It is an interesting question how the critical behaviour is modified by a confining 
geometry. If the system is finite in all directions, thermodynamic functions have 
to be analytic functions. I.e. a singular behaviour like eqs. (fT|2|) is excluded. As 
a remnant of such singularities there remains a peak in the neighbourhood of the 
transition. With increasing linear extension the hight of the peak increases and the 
temperature of the maximum approaches the critical temperature. This behaviour 
is described by the theory of finite size scaling (FSS). For reviews see |5l |6]. 

In the present work we study thin films. Thin films are finite in one direction 
and infinite in the other two directions. In this case singular behaviour is still 
possible. However the associated phase transition belongs to the two-dimensional 
universality class. I.e. in the case of U{1) symmetry, a Kosterlitz-Thouless (KT) 
transition [TJ [HI [9] is expected. In our recent work [10] we have focused on the study 
of this transition and the scaling of the transition temperature with the thickness 
of the film. 

Here we investigate the behaviour of the specific heat of thin films in the neigh- 
bourhood of the A-transition. The specific heat has been studied in a number of 
experiments on thin films of fluid ^He and ^He-^He mixtures near the A-transition. 
For recent reviews see [TTl [T2] . 

-•^In the case of the A-transition of "^He we shall denote the critical temperature by T\. 
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The physics of thin films is governed by the ratio Lq/C,, where ^ is a correlation 
length of the bulk system and Lq the thickness of the thin film. For Lq ^ ^ the 
behaviour of the film is essentially given by the thermodynamic limit of the three- 
dimensional system. In the critical region, when ^ gets close to Lq or even larger, the 
behaviour deviates from the three-dimensional one and is characterised by universal 
functions of Lq/^. 

In particular, the behaviour of the specific heat can be described by the universal 
scaling function 

[Cbuikit) - C{t, Lo)]A-'\tr = ~g2ALo/0 = g2At[Lo/^o]'^n (3) 

where Cbuik{t) is the specific heat of the three-dimensional thermodynamic limit 
and C{t, Lq) the specific heat of a film of thickness Lq. As A either ^4+ or A- can 
be taken and analogously as either ^o,+ or ^o,-- In the last step of the equation we 
have used eq. ([1]). Following RG-theory, the analytic background B is the same for 
the bulk and the thin film. Therefore it cancels in the difference that is considered 
here. Alternatively one might consider 

[CbuiM - c(t,Lo)]A-i|tr = hALo/i) = giAt[L^li^fn (4) 

where to is chosen such that ^(to) = Lq in the high temperature phase. Multiplying 
eqs. dSIlD by (|t| [Lq/^o]^/'')"" one arrives at 

nuik{t) - C{t, Lo)]A-^eo^^ = flMO = f2At[Lo/^o]'^n (5) 

and 

[C,uik{tQ) - C{t, Lo)]A-iC^^ = f\ALo/0 = fiAt[Lo/CoY^n (6) 

respectively. Often in the literature the factors A^-^ and A^-^^q^'^ are omitted and 
tL^^'^ is used as argument of the scaling function. This poses no problem as long 
as films of different thicknesses of the same system are considered. However, com- 
paring films of e.g. ^He at different pressures or ^He-^He mixtures at different 
concentrations of ^He this fact has to be taken into account. The same holds for 
the comparison of such experimental results with those obtained from lattice mod- 
els or field theory. In the following we shall use the notation /i = AC,q "^'^ fi,R and 

These universal scaling functions have been determined by a number of exper- 
iments on ^He and mixtures of ^He and ^He. In the high temperature phase, the 
data follow nicely the prediction of finite size scaling as can be seen e.g. from 
figure 14 of [12]. In this figure, [Cijuikif) — C{t, LQ)]\t\" is plotted as a function of 
[Lq/^o]^^'^ for films of ^He at vapour pressure of thicknesses 483 A up to 57/im. 
The data for different thicknesses fall nicely on top of each other. In their figure 20 
the authors of [12] have plotted data for the low temperature phase in an analogous 
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way. Up to ~ 5 the data for different tliicknesses fall nicely on top of each 

other. However for larger values of I^I^vq'''^ the data start to scatter. This is most 
pronounced at |t|i^o'^'^ ~ 10, where the function assumes a minimum. There is a 
factor of about 3.5 between the value for the thinnest and the value for the thickest 
film. From Itl-Z^^'^ ~ 20 up to ~ 100 there is a factor of about two between the 
thinnest and the thickest of the films. Note that in figure 20 of [12] Lq is given in 
A and t = 1-T/Tx. 

In the case of superfluid helium the order parameter is a complex wavefunction. 
This wave function vanishes at the boundaries of the film. In order to mimic this in 
theoretical models, Dirichlet boundary conditions with vanishing field are employed. 

Using such boundary conditions, the scaling function /2 has been calculated 
by using the e-expansion to 0(e) [131 El- The coefficients of 0(1) and 0(e) are 
numerically of similar size. Therefore one should not expect quantitatively accurate 
results for /2 obtained this way. Both fi and /2 have been computed by using 
perturbation theory in three dimensions fixed |r51|TBl|T7j in one-loop approximation. 
In the high temperature phase and at the critical point of the bulk system, the 
experimental results are fairly well reproduced. However in the low temperature 
phase, in particular close to the KT transition, no accurate predictions can be 
obtained. 

Also Monte Carlo simulations of the standard XY model on a simple cubic lattice 
have been performed to determine the specific heat of thin films. For a precise 
definition of the XY model see below. In [18] staggered boundary conditions have 
been used to obtain a vanishing order parameter at the boundaries. The authors 
of [T^ have employed free (in their notation "open") boundary conditions as we 
do in the present work. In both cases the authors have computed the finite size 
scaling function fi. In [18] the authors have simulated lattices of a thickness up to 
Lq = 24 lattice units, while in [19] the thicknesses Lq = 12, 14 and 16 have been 
studied. The results of [18] and [19] for fi agree. There is also a reasonable match 
with experiments on helium films. 

The purpose of the present paper is to compute the finite size scaling function /2 
for the first time by using Monte Carlo simulations of a lattice model. Furthermore 
we carefully study corrections to scaling, allowing us to quantify the error of our 
result for the finite size scaling function. 

For finite systems we expect that the leading corrections are cx Lq^ , irrespective 
of the type of the boundary conditions [5]. The numerical value of the correction 
exponent is = 0.785(20) for the XY universality class in three dimensions [20j : 
similar results are obtained with field-theoretic methods; see e.g. [1]. In order 
to avoid these corrections, we study an improved model. In improved models the 
amplitude of corrections oc Lq^ vanishes or in practise, it is so small that its effect 
can be ignored. The precise definition of the model that we have simulated is given 
below. 

On top of the restricted geometry, free boundary conditions introduce new phys- 
ical effects. For a discussion see e.g. reviews on surface critical phenomena [211 [22] . 
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In fact, free boundary conditions lead to additional corrections to scaling. The 
leading one is oc Lq^ [23j; it can be cast in the form Lq^j/ = + Lg. H In [10] we 
have obtained the accurate numerical estimate Lg = 1.02(7) for the model that we 
simulate here. 

This paper is organized as follows: In the next section we define the lattice model 
that we have simulated and the observables that we have computed. In section 
three we discuss how corrections caused by the free boundary conditions affect the 
finite size scaling behaviour of the specific heat. Next we discuss the details of our 
simulations. Based on these simulations we compute the scaling functions /i and 
f2- We compare our results with those from experiments on thin films of ^He, field 
theoretic methods and previous Monte Carlo simulations. 



2 The model and the observables 
2.1 The two component 0^ model 

We study the two component 0^ model on a simple cubic lattice. We label the sites 
of the lattice by x = {xq,xi,X2). The components of x might assume the values 
Xi G {1, 2, Lj}. We simulate lattices of the size Li = L2 = L and Lq <^ L. 
In 1 and 2-direction we employ periodic boundary conditions and free boundary 
conditions in 0-direction. This means that the sites with xq = 1 and xq = Lq have 
only five nearest neighbours. This type of boundary conditions could be interpreted 
as Dirichlet boundary conditions with as value of the field at Xq = and Xq = 
Lq + 1. Note that viewed this way, the thickness of the film is Lq + 1 rather than 
Lq. This provides a natural explanation of the result Lg = 1.02(7) obtained in 
[To] and might be a good starting point for a field theoretic calculation of L^. The 
Hamiltonian of the two-component (p^ model, for a vanishing external field, is given 
by 

n = -(3j2^^-^y + T.[^' + - 1)'] ' (7) 

<x,y> X 

where the field variable (px is a vector with two real components. < x,y > denotes 
a pair of nearest neighbour sites on the lattice. The partition function is given by 



n 



exp(— 7-^). 



Note that following the conventions of our previous work, e.g. [26] , we have absorbed 
the inverse temperature /3 into the Hamiltonian. In the limit A — >■ 00 the field 
variables are fixed to unit length; i.e. the XY model is recovered. For A = we get 



^ In the literature, replacing Lq by Lo^^ff — Lq + Lg to account for surface corrections, was 
first discussed by Capehart and Fisher [21] in the context of the surface susceptibility of Ising 
films. 
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the exactly solvable Gaussian model. For < A < oo the model undergoes a second 
order phase transition that belongs to the XY universality class. Numerically, 
using Monte Carlo simulations and high-temperature series expansions, it has been 
shown that there is a value A* > 0, where leading corrections to scaling vanish. 
Numerical estimates of A* given in the literature are A* = 2.10(6) [25], A* = 2.07(5) 
[26] and most recently A* = 2.15(5) [20]. The inverse of the critical temperature 
/9c has been determined accurately for several values of A using finite size scaling 
(FSS) [20]. We shall perform our simulations at A = 2.1, since for this value of A 
comprehensive Monte Carlo studies of the three-dimensional system in the low and 
the high temperature phase have been performed [IHl EOj 123, EH]- At A = 2.1 one 
gets Pc = 0.5091503(6) [20j. Since A = 2.1 is not exactly equal to A*, there are 
still corrections oc L~^, although with a small amplitude. In fact, following [20], it 
should be by at least a factor 20 smaller than for the standard XY model. 



2.2 The energy and the specific heat 

First we should note that in eq. ([7]) P does not multiply the second term. Therefore, 
strictly speaking, j3 is not the inverse temperature. However, in order to study 
universal quantities it is not crucial how the transition line in the jS-X plane is 
crossed, as long as this path is smooth and not tangent to the transition line. Here, 
following computational convenience, we vary j3 at fixed A. Correspondingly we 
define the energy density as the derivative of the reduced free energy density with 
respect to jS. Furthermore we multiply by —1 to get positive numbers: 



E 



1 d log Z 
L0L1L2 d(3 



It follows 



E 



(9) 



(10) 



We then define the specific heat as the derivative of the energy density with respect 
to /?: 

It is easy to show that 



C 



L0L1L2 



_<x,y> 



\<x,y> 




(12) 



2.3 The correlation length 

The second moment correlation length ^2nd and the transversal correlation length 
of the bulk system are used to set the scale in the high and the low temperature 
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phase, respectively. Here we shall use the results given in [28]. For completeness 
we recall the definitions of the two correlation lengths. 

The second moment correlation length in k-direction is defined by 



where 



2nd,k 



X 



xlFk - 1 
4 sin(7r/Li)2 



1 



L0L1L2 



is the magnetic susceptibihty and 

1 



L0L1L2 



E 



exp 



(13) 



(14) 



(15) 



is the Fourier transform of the correlation function at the lowest non-zero mo- 
mentum in k-direction. Note that in the high temperature phase there is little 
difference between ^2nd and the exponential correlation length ^^xp which is defined 
by the asymptotic decay of the two-point correlation function. Following [26]: 



^exp 



lim 



1.000204(3) , (t < 0) 



(16) 



for the thermodynamic limit of the three-dimensional system. |f| In [TO] we find for 
A = 2.1 by fitting the data for the second moment correlation length for f3 > 0.49: 

^2nd = 0.26362(8)(-t)°-^^^^ X [1 + 0.039(8)(-t)°-^2^ - 0.72(4)(-t)] , (17) 

where t = p- 0.5091503. 

The helicity modulus T gives the reaction of the system under a torsion. To 
define the helicity modulus we consider a system, where rotated boundary conditions 
are introduced in one direction (e.g. the 1-direction): For xi = Li and yi = 1 the 
term 4>x4'y the Hamiltonian is replaced by 

4 ■ Rjy = (cos(«)0« + sin(a)0f ) + 0(2) (_ sin(a)0« + cos(a)0f ) . (18) 

The helicity modulus is then given by 

Li dHog Z{a) 



T 



L0L2 



(19) 



Q = 



Note that we have skipped a factor of T compared with the standard definition |29j . 
Defined this way, the helicity modulus has the dimension of an inverse length. In 
the literature ^± = 1/T is referred to as transversal correlation length. Fitting the 
data for the helicity modulus at A = 2.1, given in table 2 of [28], up to /3 = 0.55, 

we get 



with t = P 



T 

0.5091503. 



1.5584(10)t' 



0.6717 



X 



;i -0.06(2)t) 



(20) 



^Throughout, in the context of the cj)'^ model we use the convention t — [3 — (3c 
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3 The finite size scaling behaviour of the specific 
heat 



In this section we discuss the finite size scahng behaviour of the specific heat of thin 
films. The free energy density of the bulk system is given by 

fbuikit) = at^-"" (1 + cit' + C2t^' + d^f + eit + . . .) + h{t) (21) 

where h{t) is the analytic background, eit is an analytic correction and Cit^, C2t^^ 
and dit^ are non-analytic corrections. In order to simplify the notation, we have 
omitted subscripts ± that indicate the phase. Numerical values of the corrections 
exponents are 9 = vuo 0.527 [20] and 9' ~ 1.2 [30]. Note that the correction 
amplitudes ci, C2, ... are small for the 0^ model at A = 2.1, while di and ei should 
assume generic values. In the following discussion we shall, for simplicity, ignore 
these corrections to scaling. Inserting ^{t) = ^ot~^ into eq. fl2Tl) we arrive at 

fbuikit) = [d^^]m-''+m (22) 

where we have used the hyperscaling relation 2 — a = du, where d is the dimension 
of the system. Note that a^g is universal. 

The free energy density of a thin film with periodic boundary conditions with 
the thickness Lq is given by [S] 

fpit, Lo) = L,Sp{t[L,/i,Y'n + m (23) 

where d is the dimension of the system and h{t) in eq. (125]) is the same function as in 
eq. ( 1211) . Also qp{x) is an analytic function at a; = 0. There might be a singularity 
at some x 7^ related with the effectively two-dimensional transition. In order to 
eliminate the analytic background one considers the difference 

huikit)-fp{t,Lo) = L^'[[d^l^]{Lo/mY-qpit[Lo/^o]'^n] 

= L^'pp{t[Lo/^ofn ■ (24) 

The specific heat is defined as minus the second derivative of / with respect to an 
analytic function of h{t). Using our definitions, it is minus the second derivative 
with respect to t itself. Applied to eq. (121]) we arrive a.t A = — /i'(0)~^(l — a)(2 — a)a 
for the amplitude of the specific heat. Taking minus the second derivative of eq. ( !24l) 
with respect to h{t) we arrive at 

auik{t)-Cp{t,Lo) = -/i(0)-%-'^[Lo/eo]'/>p(t[^o/eo]'/^) = 
-h{0)-'d[Lo/^or/''[d^X'p'Ut[Lo/^oY^n = A-'[Lo/^or^^f2,pAi[Lo/^o]'^n ■ (25) 

In the case of free boundary conditions, which are studied in the present work, 
boundary effects have to be taken into account. In [ID] we have numerically shown 
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that leading corrections can be accounted for by replacing Lq by i^o.e// = Lq + Lg 
where we have obtained the estimate = 1.02(7). Hence 



'0,eff 



(26) 



The additional term c{t)LQ^jj gives a correction of the analytic background caused 
by the free boundary conditions. Written this way it allows that the correction to 
the background is given by an effective thickness that is different from i^o,e//- The 
prefactor front of /^(t, Lq) corrects the volume that is used to compute the 

free energy density. 

Taking the same steps as in the case of periodic boundary conditions we arrive 

at 



Cbulkit) — 



(27) 



Alternatively in the literature one considers 



-Cpit, Lq) — Cbuikito) 



(28) 

o,e// for to in the high temperature phase. 



.Lo,eff 

where C(^o) = -^o or better ^(to) = L, 



Here we study the neighbourhood of the critical point. Therefore we shall ap- 
proximate w{t) ^ w{0). To simplify the notation we shall write w instead of w{0) 
in the following. 

Note that for large |t| [Lo,e///^o]^'^'' the difference Cbuik{to) - j^f^Cpit, L) be- 
comes small compared with Cbuikito) or Cp(t,L). On the other hand, corrections 
due to the boundary are virtually independent on |t| [Lo,e///Co]^^''- Therefore cor- 
rections due to the boundary might lead to large relative errors for large values of 

\t\[Lo,eff/^o]'^r 



4 Numerical Results 



4.1 Thermodynamic limit of the three-dimensional system 

In this subsection we shall consider systems with periodic boundary conditions in 
all three directions and the linear size Lq = Li = L2 = L. In the high temperature 
phase, corrections to the thermodynamic limit decay exponentially with increasing 
lattice size. In practice it turns out that for L ^ 10^2nd these corrections are 
much smaller than the statistical error that we reach here. Since a Goldstone mode 
is present in the low temperature phase of the system [311 [32] corrections to the 
thermodynamic limit decay only with some power of the linear lattice size. In 
particular for the energy density and the specific heat, we expect that, to leading 
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order, the correction decays oc L^^. Therefore rather large lattices are needed to 
get a good approximation of the thermodynamic limit. 

One should note that the estimator ( fTOl) of the energy density is self-averaging, 
while the one (|T2l) for the specific heat is not. Since the lattices have to be rather 
large to avoid sizable finite size effects, it turns out that the specific heat can be 
most efficiently determined by fitting the energy density in same range of (3. 

To this end, we have computed the energy density for a large number of (5- 
values. As starting point we have taken the data given in tables 2 and 5 of [28]. 
These were supplemented by a rather large number of new simulations to obtain 
a dense grid of /^-values. In particular, we have simulated the 96^ lattice in the 
range 0.521 < 13 < 0.58 in steps of A/3 = 0.0005, for 0.5025 < p < 0.5035 in steps 
of A/3 = 0.0001 and for 13 = 0.491, 0.492, 0.494, 0.496, 0.497, 0.498, 0.499, 0.501, 
0.502, 0.503 and 0.504. We simulated the 128^ lattice at /3 = 0.516, 0.517, 0.518, 
0.519 and 0.5205, the 192^ lattice at /3 = 0.514, 0.5145, 0.5155, 0.5165, 0.5175, 
0.5185 and 0.5195 and finally the 288 lattice at /3 = 0.5115, 0.5125, 0.5135, and 
0.5145. 

In the case of L = 96 we typically performed 10^ measurements and 5 x 10^ 
for the larger lattice sizes. For each of these measurements a Metropolis sweep, 
several overrelaxation sweeps and single cluster [33] updates were performed. For a 
discussion of the Monte Carlo algorithm see [28] . In total, these simulations took a 
little less than one year of CPU time on one core of a 2218 Opteron processor (2.60 
GHz). 

In [28j the results for the thermodynamic limit in the low temperature phase 
were obtained by fitting the data of several lattice sizes with the ansatz E{L) = 
E{oo) + cL~^. In the case of the simulations that we have added here, we have 
checked that the cL^^ corrections are sufficiently small to be ignored. 

In the neighbourhood of the transition, we have fitted the energy density with 
the ansatz 

E{(3) = Ens + Cns{(3 - f3c) + a±|/3 - /^.l^"" + d^P - f^cf + h±\f3- /^/^^ (29) 

where Ens, Cns, Pc = 0.5091503(6) and a = —0.0151(3) [20] are input and a±, dns 
and h± are the 5 free parameters of the fit. 

From the finite size scaling behaviour of systems with periodic boundary 
conditions in all directions at the critical point we find [27] : 

Ens = 0.913213(5) + 20 X (/3, - 0.5091503) + 5 x 10"^ x (1/a + 1/0.0151) (30) 

for the non-singular part of the energy density and 

Cns = 157.9(5) + 147000 x (/3, - 0.5091503) - 2.1 x {1/a + 1/0.0151) (31) 

for the non-singular part of the specific heat at A = 2.1. 

We did not include a term with the exponent l — a + uO' ^ 2 into the ansatz (!29|) . 
We expect that it is effectively taken into account by the last two terms in fl29p . 



9 



Note that the main purpose of fitting the energy density with ansatz (l29l) is to 
interpolate our data in a large range of /3-values. 

After some preliminary studies, we decided to fit the energy density in the range 
0.49 < 13 < 0.529 using the ansatz (l29l) . In total we have 98 data points in this 
interval and we get x^/d.o.f. = 1.08 for our fit. The results for the fit parameters 
are a+ = 160.688(2), a_ = -151.459(2), d^s = -302.6(9.8), &+ = 302.4(10.3) and 
6_ = 293.4(10.3). 

First we have computed the universal combination 

Ra = {l- A+/A_)/a (32) 

from the result of the fit. Note that A^/A_ = —a^/a_. Using the central values 
for the input parameters we obtain Ra = 4.035(16), where we have taken into 
account the covariance of a+ and a_. Furthermore we have checked the dependence 
of our result on the input parameters. In fact, the dependence on the value of a is 
quite small. Taking the preferred value of the experiment on the space shuttle [31] 
a = —0.0127 we get Ra = 4.022(16). We have also checked the effect of the error of 
the other input parameters. It turns out that the uncertainty of Cns has the largest 
effect on Ra. Replacing Cns = 157.9 by Cns = 158.4 we get Ra = 4.025(16). As our 
final result we quote 

i?„ = 4.035(16) [20] (33) 

where in () we give the statistical error and in [] the sum of all errors due to the 
uncertainty of the input parameters of the fit. Note that the present result is 
compatible with the final result Ra = 4.01(5) given in |27j . 

In figure [T] we have plotted the specific heat obtained from the fit of the energy 
density using ansatz fl29|) . Computing the statistical error, correlations among the 
fit-parameters are properly taken into account. In this plot, errors can not be 
resolved. In order to make the errors visible we have plotted in figure[2]the statistical 
error and the difference between the results obtained by using the central values of 
the input parameters and results where we have replaced one of the central values 
by the central value plus the error. 

As one might expect, the differences diverge in the neighbourhood of the critical 
point. In the low temperature phase the largest uncertainty is due to the error of 
Cns- In particular, going to the upper boundary (3 = 0.529 of our fit interval, the 
uncertainty induced by the error of Cns rapidly increases. Therefore we decided to 
use the values for the specific heat obtained from the fit (12^ only up to /? = 0.525. 
For larger values of (3 we follow an alternative approach as discussed below. In 
the case of the high temperature phase we shall use the results obtained from the 
fit ([29D down to /3 = 0.49. 

In order to complete the computation of the specific heat, we have fitted our 
data for the energy density with the ansatz 

n 

E{P) = Y,ai{(3-(3oy (34) 

i=0 
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0.49 0.5 0.51 0.52 



Figure 1: We plot our result for the specific heat obtained from the fit of the 
energy density using ansatz fl29|) . Note that at Pc = 0.5091503(6) the specific heat 
assumes the value Cns = 157.9(5), eq. fl?I]) . In addition to the result obtained by 
using the central values of the input parameters, we have also plotted those were 
we have replaced the central value by the central value plus the error. For example 
Cns = 157.9 is replaced by Cns = 158.4. At the resolution of the plot, all these 
curves fall on top of each other. For a more detailed discussion see the text. 
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Stat 




Figure 2: Various sources of the error of the specific heat as obtained from fits 
of the energy density using ansatz (!29l) . We plot the statistical error (stat) and 
the difference of the result using the central values of the input parameters and 
results where we have replaced the central value of one of the input parameters by 
the central value plus the error. In particular, we have replaced C^s = 157.9 by 
Cns = 158.4 (Cns), E^s = 0.913213 by E^s = 0.913218 (ens), (3c = 0.5091503 by 
Pc = 0.5091509 (be) and a = -0.0151 by a = -0.0148 (alpha). For a discussion 
see the text. 
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where ai is identified with the specific heat at /?o- We have included all values of f3 
within the interval [Pq — A, (3q + A] into the fit. We have tested various values of n. 
Our final results are taken from fits with n = 4. In the range 0.525 < jSo < 0.529 
we have used A = 0.008. For 0.529 < Po < 0.537 we have used A = 0.011 and for 
0.537 < /3o < 0.58 we have used A = 0.016. In all cases x^/d.o.f. is close to one. 
To give an impression of the accuracy that is reached, we quote C = 12.690(3) at 
p = 0.525 and C = 11.899(2) at p = 0.53. 

For comparison with the results obtained from the fit with ansatz fl29|) we have 
also performed fits with ansatz ( IMI) in the range 0.522 < /?o ^ 0.525 using A = 
0.005. The numerical results for the specific heat obtained from the two different 
approaches are compatible within error bars. 



4.2 Adjusting the scale of the axes 

As we have mentioned already in the introduction, often in the literature the fac- 
tor A~^^Q^'^ is ignored when computing the finite size scaling functions fi and f2- 
Therefore, in order to compare the results of different systems, we have to compute 

^1-2 = -— ;7- • (35) 

I ><0 \system2 

In order to fix the ratio between our definition for the specific heat of the (j)^ model 
at A = 2.1 and that of experiments on ^He at vapour pressure, we take the ratio 
of our result for A_ obtained in the previous subsection and that of [M] for the 
three-dimensional thermodynamic limit of ^He at vapor pressure. From table II of 
[51] we read off aA_^ = 5.6537. Using an alternative fit ansatz the authors get 
aA- = 5.6950. We regard this difference as an estimate of the possible error of 
From the same fits, the authors obtain a = —0.01264 and —0.01321, respectively. 
In order to match with these experimental numbers, we have taken aA_ = 2.285 
obtained from fitting the energy density with ansatz (12^ . using a = —0.0127 as 
input. Instead, using a = —0.0151 we arrive at aA_ = 2.322. I.e. the value of aA_ 
is quite insensitive on the value of a that is assumed. 

Furthermore, we need the amplitude of the correlation length for ^He at vapor 
pressure in the high and the low temperature phase. For the transversal correlation 
length in the low temperature phase one finds ^3] and refs. therein: 



^ 3.42A t-" (36) 

where t = 1 — T/T\. Alternatively we can compute the amplitudes of the correlation 
lengths from A± using the results for the universal amplitude ratios 

Rt = eo,2nd(a^+)'/' = 0.3562(10) , (37) 



*Note that our definition of A± and tliat used in [S^ differs by a factor of a 
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R- = ^o,t(«A-)^/^ = 0.850(5) (38) 

and 

= = 0.411(2) (39) 

given in [28]. To this end, we first have to convert the results for the specific heat 
of the experiment from Jmole~^if~^ into A~^. To this end we need the density 
P\ = 146.1087 kg/m^ [SB] of ^He at the A-transition and the Boltzmann constant 
h = 1.38065 ... X 10^23 g^^^j^g 

eo,T = 3.45(3)A (40) 

where we have taken into account the errors of and of the experimental estimate 
of aA-. In the case of the high temperature phase we get 

eo,2nd = 1.422(5)A (41) 

where we have used 

aA+ = 1.05251 x 5.6537 = 5.9506 (42) 
taken from table II of [31] as input and the error is estimated by using 

aA+ = 1.05490 x 5.6950 = 6.0077 (43) 



obtained from an alternative ansatz 

Now we are ready to compute the ratio 

where we can either use ^2nd or ^t- We get 1.0386 using a = —0.0151 or 1.0324 
using a = —0.0127. 
Hence we arrive at 

r.He,^, = 2.57 (45) 

with a relative uncertainty of about 2%. Let us note that this number is only 
valid for ^He at vapour pressure and the (j)'^ model at A = 2.1 and the particular 
definitions of the specific heat that have been used. 



4.3 Finite size scaling at /3c,3d 

First we performed simulations at the inverse transition temperature Pc^sd = 0.5091503(6) 
[20] of the three-dimensional system. Here the correlation length of the thin film is 
relatively small; we find C,2nd,fiim/ Lo,eff ~ 0.416. Hence already rather small ratios 
of Li/Lq are sufficient to approximate well the two-dimensional thermodynamic 
limit and therefore large values of Lq can be reached. Furthermore, by construc- 
tion Lq/C,3d = 0. Therefore this is an ideal location to accurately study the finite 
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Table 1: Results for the specific heat C obtained by using eq. / fT2j) at (3 = 0.5091503, 
which is the estimate of [20^ for (3c,3D, for lattices of the size Li = L2 = 6Lq. 





C 


8 


5.060(9) 


12 


6.145(13) 


16 


6.923(14) 


24 


8.062(17) 


32 


8.924(19) 


48 


10.141(24) 


64 


11.000(29) 



size scaling behaviour and in particular the corrections caused by free boundary 
conditions. At the critical point of the three-dimensional system eq. f l27l) reduces 
to 

-^C{t, Lo) = C^s + cLl%j + wL-l^f (46) 

where = a««fc(0) = 157.9(5) [27] and c = -^(O). 

All numbers for the specific heat discussed in this section are determined by 
using eq. ( |T2i) . As a first step, we have simulated the thickness Lq = 8 for Li = 
L2 = 16,24,32,48 and 64. We conclude from these simulations that, at the level 
of our statistical error, the two-dimensional thermodynamic limit of the specific 
heat is reached for Li = L2 = 48. Based on this result we performed simulations 
for Lq = 12,16,24,32,48 and 64 with Li = L2 = QLq throughout. In all cases, 
we performed 10^ measurements, where for each measurement we performed one 
Metropolis sweep, two overrelaxation sweeps and wall cluster |37j and single cluster 
updates. The number of single cluster updates was chosen such that the average 
size of a cluster times the number of clusters is a bit less than the number of lattice 
sites. In total these simulations took about 5 month of CPU time on one core of a 
2218 Opteron processor (2.60 GHz). The results for the specific heat, as defined by 
eq. (fT2l) . are summarized in table [1] 

We have fitted these data with ansatz (1461) . where we have fixed Cns = 157.9, 
(3c = 0.5091503, a = —0.0151 and Lg = 1.02. The results of these fits are summa- 
rized in table [2J Already for L^ rnin = 12 the x^/d.o.f. is smaller than one. Going to 
larger Lo rnin the statistical error of w rapidly increases. Therefore we take the result 
obtained for Lo^rnin = 12 as our final result. Based on our data, it is impossible to 
give an estimate of systematic errors due to sub-leading corrections. 

In order to check the dependence on the input parameters, we have repeated 
the fits for Lo,min = 12 using values of the input parameters that are shifted by 
the error of the input parameters: E.g. in one of these fits Lg = 1.02 is replaced 
by Lg = 0.95, while the other input parameters remain unchanged. In the case 
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Table 2: Fitting the data for the specific heat given in table [I] using ansatz lfW\) . 
Data for Lo^rnin < -^o < 64 are included into the fit. We have fixed Cns = 157.9, 
a = -0.0151 and = 1.02. 



f-'(i;min 


c 


w 


xVd.o.f. 


8 


-161.569(14) 


3.32(16) 


1.01 


12 


-161.597(20) 


3.90(33) 


0.21 


16 


-161.595(25) 


3.85(53) 


0.28 



of a and /3c, we have taken into account the effect of the shift on Cns and Ens as 
given by eqs. ( l30|3Ti) . It turns out that shifting Lg has the largest effect on w. 
Using Ls = 0.95 we get w = 4.40(32). Taking the shifted value Cns = 158.4 for the 
analytic background we get w = 4.17(33). Shifting the other input parameters has 
less impact on the value of w. 

Our result for c = — /2(0) can be compared with experiments and results ob- 
tained by field theoretic methods. On page 1028 in section V.A. of the authors 
analyse the scaling behaviour of the specific heat of thin films at the A-transition 
of the three-dimensional system. Fixing a = —0.01264 they arrive at 

C(0, Lo) = [453.8 ± 4.3] - [474.0 ± 4.9]Lo (47) 

where C is measured in units of Jmole'^K^^ and Lq in A. 

Analysing our data for the specific heat at i3c,3D, assuming a = —0.01264, we 
arrive at c ~ —190. Multiplying with r4He,<p-i = 2.57 we arrive at 488 in quite good 
agreement with the experimental result [474.0 ±4.9], in particular when taking into 
account the error of ri^g^^^i. 

We can also write our result in terms of the universal ratio 

- /2,ij(0) = ^ = 1.0208(2) - 1.33 X (« + 0.0151) (48) 

where the error quoted in () is dominated by the statistical error of the specific heat 
of the thin films and the uncertainty of Ls- The dependence on the value of a that 
is used in the analysis is rather weak. The authors of [13] have computed — /2,_r(0) 
(in their notation uJoo) using e-expansion to 0(e). Their result is given in their 
eq. (8.17). To leading order ujoo = 1, which is in quite good agreement with both 
our numerical result and with the experiment. Dohm and Sutter [T7] have pointed 
out that the extrapolation to e = 1 is affected by considerable ambiguities. In fact, 
setting e = 1 the authors of [13] find uoo ~ 0.75, (table II of [13]) which is clearly 
ruled out by our result as well as by experiment. 
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Table 3: Results for the speciGc heat at the KT transition of thin films. In the 
first column we give the thickness Lq of the film, in the second column, we give the 
result for Pkt{Lo) obtained in In the third column we give the estimate of the 
specific heat Cbuik of the three-dimensional bulk system at this temperature and 
finally, in the fourth column the specific heat of the thin film at the KT transition. 
Throughout we have taken the result for Li = L2 = 32Lo. For a detailed discussion 
see the text. 





I3kt 


Cbulk 


C film 


6 


0.56825(1)[1] 


8.586(1) 


9.710(15) 


8 


0.549278(5) [9] 


9.870(1) 


10.793(17) 


12 


0.532082(3) [5] 


11.617(2) 


12.179(20) 


16 


0.524450(2) [3] 


12.791(7) 


13.176(23) 


24 


0.517730(2) [2] 


14.342(5) 


14.454(22) 


32 


0.514810(1)[2] 


15.379(6) 


15.311(33) 



4.4 Finite size scaling at I3kt{Lq) 

The KT phase transition occurs, up to scahng corrections, at a given value of the 
scaling variable thU" or equivalently Lo/C,t- In [lO] we find [Lq/^t]* = 1.595(7). 
Following the KT theory [71 El IS] the free energy is infinitely often differentiable with 
respect to the temperature at the KT transition. Therefore no particular problem 
for the numerical analysis is expected. At the KT transition eq. fl27j) reduces to 

CbuikiPMLo)) - -^CiPMLo), Lo) = CKT^Teff - wL,l.. (49) 
Wff 

where ckt is /2 at [Lq/^t]* = 1.595(7). 

We have obtained accurate data for several values of Lq in relation with [10]. 
In table [3] we have summarized the results for the specific heat obtained by using 
eq. (fT2|) at the KT transition. The results that are quoted were obtained for Li = 
L2 = 32Lo lattices. In the case of Lq = 24 and 32, these are the largest available. 
For Lq < 24 we have checked that the results obtained from Li = L2 = 32Lo 
are consistent within the statistical error with those obtained from larger values 
of Li = L2. In addition, we give numerical estimates for the specific heat of the 
three-dimensional bulk system as discussed in section 14. 1[ 

We have fitted these data with ansatz (H9|) . Our results, using the central values 
of the input parameters, are summarized in table HI The value of w is increasing 
with increasing minimal thickness Lo^rnin that is included into the fit. The result 
obtained for Lo^min = 16 is compatible within the statistical error with that obtained 
in the preceding subsection for Pc,3D- 

Next, we have checked how the results for ckt and w depend on the values of the 
input for Lg and a. To this end we have repeated fits for Lo^min = 16 using shifted 
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Table 4: Fits of the specific lieat C at tie KT transition witli ansatz ( BPj) . Data 
with Lo^rnin < -^0 < 32 are included into the fit. We have fixed a = —0.0151 and 
Ls = 1.02. 



-f-'0,min 


CRT 


w 


xVd.o.f. 


8 


0.587(21) 


2.07(19) 


2.98 


12 


0.638(27) 


2.81(31) 


1.00 


16 


0.650(42) 


3.04(66) 


1.42 



values for these input parameters. Using a = —0.0151 and Ls = 0.95 as input we 
get ckt = 0.634(42) and w = 3.61(65). Changing the value of the exponent of the 
specific heat to a = —0.0148 and keeping Lg = 1.02 we get ckt = 0.649(42) and 
w = 3.03(66). Finally we have set a = —0.0127 obtained from the space shuttle 
experiment [M]. Keeping Lg = 1.02 we get ckt = 0.641(41) and w = 2.99(65). 

As one might expect, w shows a strong dependence on the value of Lg. On the 
other hand the dependence of w on a is quite weak; even using the preferred value 
of [31] the results for ckt and w change only little. As our final estimate we take 
Ckt = 0.65(4) [2]. Taking into account both the analysis at Pc,3D and at Pkt we 
shall use w = 3.5 in the following analysis of the specific heat in a large range of 
the scaling variable. The error of w should be about 1. 

4.5 Specific heat for Lq = 8, 16 and 32 for a large range of /3 

Finally we have computed the specific heat for Lq = 8, 16 and 32 for a large range of 
P in the neighbourhood of Pc,3D- For this purpose, it turns out to be more efficient 
to compute the specific heat by taking the derivative of the energy density (fTOll 
with respect to (3 numerically than by using eq. (fT2l) . 

Let us first discuss how we have computed the two-dimensional thermodynamic 
limit of the energy density of the thin films. In the high temperature phase of 
the thin film we expect that the energy density converges exponentially fast with 
increasing L = Li = L2 toward the effectively two-dimensional thermodynamic 
limit. In order to see this asymptotic behaviour, lattices with L ^ C,2ndjum are 
needed. From numerical results for the two-dimensional XY model we conclude 
that L ^ 8^2nd,fiim is needed such that the deviation from the thermodynamic limit 
is by far smaller than the statistical error that we typically reach in our study. 
For Lq = 8, 16 and 32 we have simulated lattices up to L = 2048, 1800 and 
1024, respectively. Therefore we could satisfy the condition L ^ S^2ndjiim up to 
P = 0.545, P = 0.522 and (3 = 0.5134, respectively. At these values of [3 we find 
^2ndjiim = 242.5(2), 153.16(15) and 111.99(15) for films of the thickness Lq = 8, 16 
and 32, respectively. 

For j3 > (3kt the asymptotic behaviour is given by the spin-wave approximation; 
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i.e. by a free field theory. Tlierefore the thermodynamic hmit is approached as 
E{L) = E{oo)+0{L-^). In this range of/5 we have taken E(oo) = |[4E(2L)-E(L)] 
as our final result for the thermodynamic limit. 

Unfortunately there is a quite large range of (3, where the extrapolation to the 
thermodynamic limit is less clear: Pmax < (3 < Pkt, where (3max is the largest value 
of P such that L^ax ~ ^^2nd and L^ax is the largest lattice size that we can simulate 
with our (finite) computer resources. 

In the context of [38] and ^9] we have simulated the two-dimensional XY model 
at /3 = 1.1199, which is the best estimate of the inverse KT transition temperature 
|40] . on lattices up to L = 4096. Fitting the energy density for 16 < L < 4096 with 
the ansatz 

E{L) = E{oo) + cL-' (50) 

where E{oo), c and e are free parameters, we get e = 1.85(2) with x^/d.o.f. smaller 
than one. Here we do not intend to further discuss this phenomenological observa- 
tion; I.e. whether this effective exponent is e.g. caused by logarithmic corrections. 
We have also generated data for the 2D XY model for various values of /3 in the 
range Pmax < 13 < Pkt for various lattice sizes L with L ^ 4^2nd- These data 
can be nicely fitted with the ansatz ( l50l) . where now e apparently depends on /3; 
it is decreasing with decreasing /?. E.g. for /? = 1.0929 we find, fitting the data 
for 16 < L < 512, the effective exponent e = 1.67(2). And for f3 = 1.0, where 
^2nd = 40.09(8), we find, fitting the data for 16 < L < 128, the effective exponent 
e = 0.93(1). 

From these observations in the two-dimensional XY model we learn that for 
Pmax < /3 < Pkt the extrapolation to the thermodynamic limit is non-trivial. In 
order to keep the systematic error small, we simulated as large lattices as possible. 
Since we had only few lattice sizes at hand, and the effective exponent e is not 
known a priori, we have also extrapolated our data in the range Pmax < /3 < Pkt 
using e = 2. In order to estimate the systematic error of this extrapolation, we 
compared results obtained for different L. 

In table[5]we have summarized the lattice sizes Lq, L = Li = L2 and the values of 
(3 that we have simulated. Typically we have performed 10^ measurements for each 
simulation. Details are given in table O For each measurement we performed one 
Metropolis sweep, several overrelaxation sweeps, single and wall-cluster updates. 
Integrated autocorrelation times in units of measurements for the energy density 
are TE,int ~ 10 for these simulations. Note that in the case of Lq = 32 we have 
skipped the range Pmax ~ ~ Pkt, since we were not able to simulate sufficiently 
large L to get good control on the thermodynamic limit. In total these simulations 
took about three years of CPU-time on a single core of a 2218 Opteron processor 
(2.60 GHz). 

Similar to the case of the three-dimensional bulk system in the low temperature 
phase, we have computed the specific heat by fitting the data of the energy density 
with the ansatz (IMIl . Also here we have taken our final results from fits with n = 4. 
We have adjusted the range of the fit A such that x^/d-o.f. is about one. 
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Table 5: We give a compilation of the simulations that we have performed. In the 
first column we give the thickness Lq of the film. In the second column we give 
the size L = Li = L2 of the lattice in the other two directions. In the third and 
fourth column we give the upper and lower bound of the interval in the inverse 
temperature that has been simulated. In the fifth column we give the step size A/3 
that has been used. E.g. /3^j„ = 0.49, P^ax = 0.52 and A/5 = 0.001 means that 
(3 = 0.49, 0.491, 0.492, . . ., 0.52 have been simulated. Finally, in the last column we 
give the number of measurements (stat) that have been performed for each of the 
simulations. 



T 


T T 

Ijl — 1^2 


Q 

Pmin 




Pmax 


A/? 


stat 


8 


64 


0.49 


0.52 


0.001 


5 X 10^ 


8 


128 


0.52 


0.527 


0.001 


2 X 10^ 


8 


128 


0.544 


0.56 


0.001 


2 X 10^ 


8 


256 


0.528 


0.537 


0.001 


10^ 


8 


256 


0.544 


0.56 


0.001 


10^ 


8 


256 


0.562 


0.58 


0.002 


10^ 


8 


512 


0.5375 


0.548 


0.0005 


10^ 


Q 


oiz 


U.o4y 


U.OD 


U.UUi 


1 n5 
iU 


8 


512 


0.562 


0.58 


0.002 


10^ 


8 


1024 


0.5435 


0.548 


0.0005 


10^ 


8 


2048 


0.545 






14 X 10^ 


16 


128 


0.49 


0.4995 


0.0005 


10^ 


16 


256 


0.527 


0.550 


0.001 


2 X 10^ 


16 


512 


0.500 


0.512 


0.001 


10^ 


16 


512 


0.5125 


0.529 


0.0005 


10^ 


16 


512 


0.53 


0.55 


0.001 


10^ 


16 


1024 


0.521 


0.529 


0.0005 


10^ 


16 


1600 


0.5223 






10^ 


16 


1800 


0.522 


0.523 


0.00025 


5 X 10^ 


32 


256 


0.49 


0.505 


0.0005 


10^ 


32 


256 


0.50525 


0.50875 


0.00025 


10^ 


32 


256 


0.516 


0.518 


0.0002 


10^ 


32 


256 


0.5185 


0.526 


0.0005 


10^ 


32 


512 


0.509 


0.5128 


0.0002 


2 X 10^ 


32 


512 


0.516 


0.518 


0.0002 


10^ 


32 


512 


0.5185 


0.526 


0.0005 


10^ 


32 


1024 


0.513 


0.5136 


0.0002 


5 X 10^ 
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In the case of Lq = 8 we have used in the neighbourhood of the maximum of 
the specific heat A = 0.003. As we go away from the maximum A is increased up 
to A = 0.011 for the smallest and largest values of f3 that we have simulated. 

In figure [3] we plot results for the specific heat for Lq = 8 in the most difficult 
range of (3 around the maximum of the specific heat. With this plot we like to check 
two sources of systematic error: First the truncation effects of eq. (1341) . To this end 
we have plotted the results obtained with = 3 in addition to those obtained with 
n = 4. Second we investigate systematic errors of the extrapolation to the two- 
dimensional thermodynamic limit in the range f3max < < Pkt- To this end we 
have replaced the values for the energy density obtained from the extrapolation of 
L = 512 and L = 1024 (set 1) by those obtained from the extrapolation of L = 256 
and L = 512 (set 2). 

The results obtained from n = 3 and n = 4 fall nicely on top of each other. 
Hence there should be no systematic errors due to the truncation of eq. (I34p that 
are considerably larger than the statistical error. On the other hand, comparing the 
results obtained from set 1 and set 2 we see discrepancies that are a few times the 
statistical error. Note that this problem affects only the range Pmax < /5 < Pkt- 
The statistical error of the specific heat is maximal at the peak of the specific heat. 
There it is about 0.01. Instead, computing the specific heat by using eq. ( fT2|) we 
get e.g. for P = 0.543, L = 512 the result C = 12.52(9), for /? = 0.5435, L = 512 
the result C = 12.62(9) and for /? = 0.5435, L = 1024 the result C = 12.48(11). 
I.e. the values for the specific heat are compatible with those obtained by fitting 
the energy density, however the statistical error is about ten times larger. 

In the case of Lq = 16 we have used in the neighbourhood of the peak of the 
specific heat A = 0.001 for the fits of the energy density. As we go away from 
the maximum A is increased up to A = 0.015 for the smallest and largest values 
of (3 that we have simulated. At the maximum of the specific heat the statistical 
error is about 0.01. The ratio L/Lq that we have maximally reached for Lq = 16 
is smaller than for Lq = 8. It is comparable with that of the data included into 
set 2 for Lq = 8 discussed above. Therefore we expect that deviations from the 
two-dimensional thermodynamic limit in the range (3max < ~ Pkt are of similar 
size as for set 2, i.e. a few times the statistical error that we have reached. Note 
again that outside of this interval, the two-dimensional thermodynamic limit is well 
under control. 

In the case of Lq = 32 we did not study the range Pmax < ~ Pkt, since 
we could not simulate sufficiently large lattices to get a good approximation of the 
two-dimensional thermodynamic limit. For the /5-values closest to the peak of the 
specific heat, we have used A = 0.0015. For our largest and smallest values of /3 we 
have used A = 0.015. Close to the peak the statistical error of C is about 0.01. 

In figure H] we have plotted our results for the specific heat for the three- 
dimensional bulk system and the thicknesses Lq = 8, 16 and 32. 

The position of the peak of the specific heat approaches the transition tempera- 
ture of the three-dimensional bulk system and the hight of the peak increases as the 
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Figure 3: In the figure we give results for tlie tliermodynamic limit of the specific 
heat defined by eq. (fTTj) for films of the thickness Lq = 8. These results were 
obtained by fitting our data for the energy density with the ansatz using n = 3 
and n = 4. We have used two different sets (set 1, set 2) of data for the energy 
density. For a detailed discussion see the text. Note that Pkt = 0.549278(5) [9] for 
Lq = 8, while the maximum of the specific heat is located at /3 ~ 0.5432 as we read 
off from the plot. 
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Figure 4: We plot our results for the specific heat as a function of the inverse 
temperature (3. In addition to the 3D bulk specific heat we give our results for the 
thicknesses Lq = 8, 16 and 32. 
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thickness Lq of the film increases. It is interesting to note that in the case of Lq = 8 
for low temperatures the specific heat of the film is larger than for the bulk system, 
while for Lq = 32 it is smaller. The specific heat of the films at the maximum is 
larger than the bulk specific heat. 

5 Results for the finite size scaling function /2 

In this section we compute the finite size scaling functions /2 using the numerical 
data for the specific heat discussed above. Our results are compared with those 
obtained from experiments on films of ^He near the A-transition and with results 
obtained by using field theoretic methods. 

5.1 The high temperature phase 

First we compute /2 without taking into account corrections. To this end, in figure 
[5l we have plotted [Cbuikif^) — C{l3, Lq)]Lq°'^'^ as a function of Lo/^2nd(/?), where 
C2nd{P) is given by eq. (fT7|) . Corrections are clearly visible. Note that in the whole 
range of Lo/^2nd{P) that is plotted, the error of [CbuikiP) — C{P,Lq)] should be at 
most 0.01. I.e. the error is much smaller than the difference between the results 
obtained for different Lq. 

Next, in figure [6], we have taken into account boundary corrections by replacing 
Lq by Lo^eff = Lq + Lg with Lg = 1.02. This has two effects: On the x-axis we 
replace Lo/^2nd by i^o,e///^2nd and secondly, in order to compute the energy density, 
we replace the volume LqLiL2 by the effective volume Lo,e//-^^i-^^2- Hence the specific 
heat that we have computed before is multiplied by Lo/Lo,e//- Now the corrections 
are much reduced in comparison with figure [51 In particular the curves for Lq = 16 
and Lq = 32 fall almost on top of each other. 

Finally, in figure [71 we have taken into account corrections due to possible 
boundary effects on the analytic background by adding the term w/Lo,e//- The 
curves match only marginally better than in figure [6l Note that while in the 
case of figure [6] the value of /2 was increasing with increasing Lq, it is decreas- 
ing in the case of figure U\ Therefore it seems to be reasonable to assume that the 
asymptotic result for /2 is located between [CbuikiP) ~ -^o)-^o/-^o,e//]-^oe// and 

[Cbuik{f3) - C(/3,Lo)Lo/Lo,e// + w/Lo,e//]^(^e/J with w = 3.5. 

In table [6] we give /2 for a few values of Lq^^/j /(,2nd- This should help the reader 
to compare our result with that obtained from other systems. 

In figure [H] we plot [Cbuik{t) — C{t, Lq)]Lq°'^'^ as a function of LQ/^2nd{t) using 
the experimental data of [HI [12] for thin films of ^He at vapor pressure of the 
thicknesses 483, 1074, 2113, 5039, 6918 and 9869A. These data are taken from the 
web page |13]. The experimental data for the specific heat are given as a function 
of the reduced temperature t = (T — T\)/T\. In order to plot them as a function of 
LQ/C2nd(t) we use (2nd = 1.422 jtl'^^A, eq. ( 1411) . As value of the critical exponents 
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Figure 5: We plot [ChuikiP) — C{P, Lq)]Lq"''^ as a function of LQ/^2nd{P) for the 
thicknesses Lq = 8, 16 and 32. For a discussion see the text. 



Table 6: We give our results for /2 obtained from Lq = 32 and corrections charac- 
terized by Lg = 1.02 and w = 3.5 for a few values of LQ eff/^2nd (upper row). Using 
instead Lg = 1.02 and w = 0, the value of f2 is smaller by about 0.11 throughout. 
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Figure6: We plot [Cbuik{P)-C{P, Lo)Lo/ Lo^effjL^^ff as a function of Lo,e///6nd(/5) 
for the thicknesses Lq = 8, 16 and 32. For a discussion see the text. 



26 




Q I I I I I I I I I I I I 

1 2 3 4 5 6 

Lo,eff^^2nd 

Figure 7: In the figure we plot [Cbuik{/3) - C{f3, Lo)Lo/Lo,e// + w/-^vo,e//]-^^(^e// with 
w = 3.5 as a function of i^o,e///^2nd(/3) for the thicknesses Lq = 8, 16 and 32. For a 
discussion see the text. 
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Figure 8: Experimental results 1^ Hlj obtained for films of ^He in the high 
temperature phase, close to the A-transition. We plot [Cbuik{t) — C{t, Lq)]Lq°'^'^ as a 
function of Lo/6nd(t) for the thicknesses Lq = 483, 1074, 2113, 5039, 6918, 9869A 
and 57 /im. For comparison we give our results obtained for Lq = 32. In contrast 
to the experimental results, boundary corrections, characterized by = 1.02 and 
w = 3.5 (solid black line) or by = 1.02 and w = (dashed black line) are taken 
into account. Furthermore we have multiplied our results for /2 by r4//e,</>4 = 2.57. 
For a discussion see the text. 

we take u = 0.6717 [20j and correspondingly a = —0.0151. Furthermore we plot 
the results of [S] for a film of 57 fim as thickness. The numbers [ISj used for the 
plot are those of [H] plotted in figure 29. We have replaced x = tLH" by i^o/^2nd 
on the X-axis and have multiplied /a by a factor ^^{){)(f-m5i/omn-Q.m27/Qmm = 104. 
Note that the authors of fTl] assume v = 0.6709 [33], while we prefer u = 0.6717 

m- 

For comparison we give our result for /a obtained from Lq = 32, taking into 
account boundary corrections characterized by Lg = 1.02 and w = 3.5 (solid black 
line). To indicate the possible error we give in addition the result obtained for 
Lq = 32, Lg = 1.02 and w = (dashed black line). We have multiplied our 
numbers by r4H^ ,f,4 = 2.57, eq. (1451) . 
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We observe that /2 computed from the specific heat of 483, 1074 and 2113 A 
films is systematically larger than our result. In contrast, we see a quite good match 
with the results obtained from 5039 and 6918 A and (a little worse) for 9869A films. 
Also in the case of the 57 /im film we see a reasonable match with our result. Note 
that there are experimental results available for much larger and for smaller values 
of Lq/C, than plotted in figure [81 

The scaling function /2 has been calculated perturbatively to one-loop in three 
dimensions fixed [151 [16]. Also in this case, the relative factor between the experi- 
mental and theoretical results for the specific heat was fixed by using the behaviour 
of the specific heat in the thermodynamic limit. To this end the experimental data 
of [l6] had been used. In [15] the result for /2 is only given as log- log plot. There- 
fore we abstain from plotting it in figure [H E.g. in figures 29 and 31 of [IT] the 
field theoretic result of [15] is plotted along with the experimental results of |34j . 
There is a quite reasonable match between the field theoretic result and the exper- 
imental data. For tL^^'^ ^ 12, the experimental result for /2 is somewhat larger 
than that of To pick out one point: For tLQ^'^ = 1 we read off from figure 

31 of [llj for [15] the value /2 ~ 10. This has to be compared with our result 
f2{i2nd/Lo = 1/1.422) = 12.6(1), where the factor riHe,(P'i = 2.57 has been taken 
into account. 

In eq. (8.1) of [13] the specific heat of the thin film has been calculated up to 
0(e) as a function of Lo/^2nd- Here we abstain from evaluating this function. In 
section [T3l we have discussed the case Lo/^2nd = 0. Below, in section [531 we shall 
discuss the limits Lo/^2nd oo and Lq/^t oo. 

5.2 The low temperature phase 

Next we did the same exercise using our data in the low temperature phase. First we 
compute the finite size scaling function /2 without taking into account corrections. 
To this end, in figure [H we have plotted [Cbuikif^) — C{(3, Lq)]Lq°'^'^ as a function of 
Lq/^t{P), where T(/3) = 1/^t{P) is given by eq. (12U1) . As in the high temperature 
phase, scaling corrections are clearly visible. Deep in the low temperature phase, 
the function even changes the sign as the thickness Lq of the film increases. Note 
that these differences are clearly larger than the statistical errors. As we have 
discussed before, the error of [CbuikiP) — C{i3,Lq)] should be at most 0.01 outside 
of the interval Pmax < P < Pkt and maybe up to 0.05 inside of this interval. This 
almost directly translates into the error of [Cbuikif^) — C{(3, Lo)]L'^°'^'^ since L^"^'^ 
is close to one for Lq = 8, 16 and 32. 

In figure [TOl we have taken into account the leading corrections to the singular 
part of the specific heat by replacing Lq by -Z^o,e// = Lq + Lg. This has two effects: 
On the X-axis we replace Lq/Ct by Lo,e///^r and the specific heat of the film is 
multiplied by Lo/-Lo,e//- After this replacement, the three curves are much closer 
than in figure[9l Note that here the error induced by the uncertainty of Lg = 1.02(7) 
dominates the error of [CbuikiP) - C(/5, Lo)Lo/-Lo,e//]-^^(^e//- 
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Figure 9: We plot [Cbuik{(3) — C(/3, Lo)]-^o "^'^ as a function of Lo/C,t{P) for the 
thicknesses Lq = 8, 16 and 32. For a discussion see the text. 
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0,ef/^T 



Figure 10: We plot [Cbuik{P)-C{P, Lo)Lo/ Lo^eff]LQ^ff as a function of Lo,e///^T(/5) 
for the thicknesses Lq = 8, 16 and 32. For a discussion see the text. 
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Figure 11: We plot [Cbuik{f3) - C(/3, Lo)Lo/Lo,e// + w/Lo,e//]-^^(^e/j with w = 3.5 
as a function of Lo^eff /^t{P) for the thicknesses Lq = 8, 16 and 32. For a discussion 
see the text. 

Finally, in figure [TT] we have taken into account the correction w/Lo,e//, where 
we have set w = 3.5. Now the curves fall nicely on top of each other giving 
support to the suggestion made in section [3] that boundary correction to the analytic 
background are not well described by replacing Lq by Lo^eff = Lq + Lg. Note that 
now the uncertainty of w gives the largest contribution to the error of [CbuikiP) — 

In tables [7] and [8] we give /2 for a few values of LQ^ff/^T- This should help the 
reader to compare our result with that obtained from other systems. 

The finite size scaling function shows a clear minimum at a finite value of 
[Lo,e///^T]- By construction, the position of this minimum does not depend on 
w. For Lq = 16 we get [i^o,e///^T]mm ~ 1.452. For w = 3.5 the minimum takes the 
value —0.47 and for = the value —0.69. For Lq = 8 we get [-C/o,e///Cr]min ~ 1.468 
and as value —0.36 for w = 3.5 and —0.76 for w = 0.0. We consider it as a robust 
result that /2 assumes a negative value at its minimum. Note that the KT transi- 
tion takes place at Lq/C,t = 1.595(7) [TU]. I.e. the minimum of /2 is located at a 
temperature slightly higher than the transition temperature. 
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Table 7: We give /2 for a few values of LQ^eff/ir (upper row). These numbers are 
obtained from Lq = 32 and corrections characterized by Ls = 1.02 and w = 3.5. 
Assuming Ls = 1.02 and w = 0, the values of f2 are lower by about 0.11 throughout. 
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Table 8: We give f2 for a few values of Lo^eff/^T (upper row). These results are 
obtained from Lq = 16 and corrections characterized by Lg = 1.02 and w = 3.5. 
Assuming Lg = 1.02 andw = 0, the values of /2 are lower by about 0.22 throughout. 
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In figure [l2] we compare our result for /2 with experimental ones. To this end, 
we plot [Cbuikit) — C(t, Lq)]Lq°'^'^ as a function of Lo/^t(^) using the experimental 
data of [m |l2] for thin films of ^He at vapor pressure of the thicknesses 483, 
1074, 2113, 5039, 6918 and 9869A. These data are taken from the web page [43]. 
The experimental data for the specific heat are given as a function of the reduced 
temperature t = (T — Tx)/Tx. In order to plot them as a function of Lo/^Tif) 
we use = 3.45 jtl^'^A, eq. (HOj) . As value of the critical exponents we take 
V = 0.6717 [20] and correspondingly a = —0.0151. Furthermore we plot the results 
of [H] for a film of 57 fim as thickness. We have rescaled these data as discussed 
in the previous subsection on the high temperature phase. For comparison we give 
our result for /2 obtained from Lq = 16 and 32, taking into account boundary 
corrections characterized by Ls = 1.02 and w = 3.5 (solid black lines). To indicate 
the possible errors we give in addition the results obtained for Ls = 1.02 and w = 
(dashed black lines). We have multiplied our numbers by Vi^^^^i = 2.57, eq. fj^5|) . 

Here almost all experimental results for the finite size scaling function /2 are 
somewhat larger than ours. There is some scattering among the experimental re- 
sults. Those for the larger thicknesses of the film are closest to our /2. There is 
a nice match for the position of the minimum of /2 obtained from 57 fim and our 
result. However /2, computed by using the experimental data, never assumes neg- 
ative values. Also the dip around the minimum is much less pronounced than it is 
in our case. It is beyond the scope of the present work to discuss possible sources 
of these discrepancies. This requires detailed discussions with the experimentalists. 
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Figure 12: Experimental results [ITl US, HI] obtained for films of ^He in the low 
temperature phase, close to the A-transition. We plot [Cbuik{t) — C{t, Lo)]Lq'^^'^ 
as a function of Lo/^rit) for the thicknesses Lq = 483, 1074, 2113, 5039, 6918, 
9869A and 57 /xm. For comparison we give our results obtained for Lq = 16 and 
32. In contrast to the experimental results, boundary corrections, characterized by 
Lg = 1.02 and w = 3.5 (solid black lines) or by Lg = 1.02 and w = (dashed black 
lines) are taken into account. Our results for /2 are multiplied by r4He^^4 = 2.57. 
For a discussion see the text. 
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5.3 The finite size scaling function /2 in the neighbourhood 
of the critical point 

Here we like to give an explicit formula for /2 in the neighbourhood of t = 0. To 
this end, we use the results for the amplitudes A± obtained in subsection I4.H the 
amplitudes of the correlation length ^o,2nd and ^o,t and the value of /2(0) obtained 
in subsection 14.31 In addition we have extracted from our data the slope of the 
specific heat of thin films at the transition temperature of the three-dimensional 
bulk system. 

As result we get in the high temperature phase: 

h = 161.6(6) -2.1(5 -[158.3(6)-25](Lo/6nd)""/'^ + 0.65(2)(Lo/6nd)'/'^ + ... (51) 

where 6 = l/a + 1/0.0151 gives the dependence on the value of a that is used for the 
analysis. This formula gives a good approximation of /2 up to about i^o/^2nd = 0.7, 
where the deviation from the full result as computed above is less than 1%. 
In the low temperature phase we obtain: 

h = 161.6(6) - 2.15 - [152.2(6) - 25](Lo/eT)""/'^ - 2.43(5) (Lq/Ct)'/'^ + • • • (52) 

This is a good approximation of /2 up to about Lq/^t = 0.3, where the deviation 
from the full result as computed above is less than 1%. 

5.4 The finite size scaling function /2 for large \t\ 

In the high temperature phase, for Lq ^ ^2nd,3D, the two boundaries are uncorre- 
lated and therefore the dependence of physical quantities on Lq is trivial. In the 
case of the specific heat we can write (ignoring boundary corrections) 

LoC(/3, Lo) = LoCtuikW) + 2Cs{f3) for (3 < (3, , Lq > 6nd,3D • (53) 

Inserting this equation in the definition of the scaling function /2 we get 

/2 = -2C,(/3)Lo for Lq > 6nd,3D (54) 

Since /2 is a function of Lo/^2nd only and Cs{l3) does not depend on Lq, it follows 

Cs = h^e,:^"' (55) 

following the convention in the literature 

= ' ^"'^"'^ (56) 
a + u 

In the low temperature phase, the situation is more complicated, due to the presence 
of the Goldstone mode. Power like corrections are present in this case: 

LoC(/3, Lo) = Lo[Cbuik{f3) + 0(1^')] + 2C,(/3) for (3 > (3, , Lq > ^t.sd • (57) 
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We do not have data for a sufficiently large range of temperatures to check carefully 
this behaviour. Assuming the correctness of eq. (15^ we read off from our data for 
Lo = 32: 

h+ = -2.2(2) 



in the high temperature phase and 

/i_ = -1.2(2) 
in the low temperature phase, where 



h. 



a + p 



(58) 
(59) 

(60) 



In both case, we have taken into account the boundary corrections. It follows that 
= —5.2(5) and A~ = —6.8(1.1) in units of jAmole~^i^'~^. Furthermore, we can 
compute the universal ratio 



Q = ^ = 0.8(2) . 



(61) 



Experimental results for films of He have been summarized by the authors of [12 
as 



At 



-5.9 ±0.2 



A: 



-8.6 ±0.5 



in units of jAmole ^. It follows 

Q 



= 0.69(5) 

At 



(62) 



(63) 



Our result for A'^ is consistent with that of experiments. In the case of A~ we see 
a small discrepancy. The results for Q are consistent within the error bars. 

The surface specific heat has been calculated using the perturbative expansion in 
three dimensions fixed in the two- loop approximation by Mohr and Dohm |17l HE]- 
Inserting numerical values into eqs. (3,4) one gets A'^ = —5.429 and Aj = —1.822 
in units of JAmo\e~^ [19]. We notice that the value for Af is in excellent 
agreement with our results and also close to the experimental one. In contrast, 
the value for Aj is clearly ruled out by us as well as by the experiment. Already 
the authors of [57] have pointed out that their result for Aj does not provide an 
accurate numerical estimate. 

Krech and Dietrich quote as result of the e-expansion, eq. (E6) of 



A + td-l 



2567r 



2±e 



2 ± In vr - 7 ± 



A^±2' 

ivTs 



+ 0(6^) 



(64) 



Inserting N = 2 one gets 



-0.0049736 - 0.0073796e± 



(65) 
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In order to compare with the experiments on superfluid '^He at vapour pressure we 
insert = 1A22A and then convert from units of to jA/mole~^i^~^. We 
arrive at 

A+ = [-1.88 - 2.79e + . . .] jA/mole"^ir-^ (66) 
Following Eisenriegler [50j 



A+ AT 

where we have inserted N = 2. The result of Eisenriegler has been extended by 
Upton [51] to 0(e2). For = 2 one gets: 

= 0.222e[l + 0.93e + . . .] (68) 
For a detailed discussion of the field theoretic results see the PhD thesis of Mohr 



6 The finite size scaling function fi 

In order to compute the finite size scaling function fi we have to determine the 
specific heat at the temperatures, where ^2nd of the three-dimensional bulk system 
assumes the values 8, 16 and 32 or, taking into account corrections, 9.02, 17.02 and 
33.02. To this end, we have first numerically inverted eq. ( fT7|) . The corresponding 
values of /? are 0.50295, 0.50694, 0.50836 and 0.50396, 0.50713, 0.50840, respectively. 
At these values of (3 the specific heat of the three-dimensional bulk system assumes 
the values 7.091(2), 9.277(2), 11.512(4) and 7.462(2), 9.475(2), 11.613(4), respec- 
tively. These values are obtained by using the results of fits with the ansatz (!29|) . 
The error is dominated by the uncertainty of Cns that has been used as input in 
eq. m- 

In figure [12] we plot [j^^^C{t, Lq) — Cbuik{to) ~ '^/-^o,e//]-^oe// ^ function 
of t{LQ^eff/io,2ndY^''i whcrc wc havc used Lg = 1.02 and w = 3.5. In particular 
for the low temperature phase we see a strong dependence on Lq. We argue that 
this is due to analytic corrections. In the case of /2 these corrections affect the 
specific heat of the three-dimensional bulk system and the thin film in the same 
way. Therefore there is only a rather small effect on the difference of the two. In the 
case of /i the specific heat of the three-dimensional bulk system and the thin film 
are taken at different reduced temperatures. Therefore there is a quite huge effect 
on the difference of the two. Motivated by this, in figure [HI we have subtracted 
Cbuikito) -|- 30 X {Po — (3) instead of Cbuik(to). The coefficient 30 is fitted by eye to 
get a reasonable collapse of the curves at low temperatures. 

In figure [15] we compare with experimental results for the finite size scaling fi. 
In order to compute to we have used ^2nd = l-4:22At~^-^'^^'^ , eq. (HTj) . We have plotted 
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t(Lo,eff/^o) 



Figure 13: We plot [j^C{t,Lo) - C^uikito) - w/ Lo^eff]L(^^'/f versus 

t(Lo,e///^o,2nd)^^'^, where we use Lg = 1.02 and w = 3.5. For a discussion see 
the text. 
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Figure 15: In the figure we plot [C(t, Lq) — C;,„;fc(^o)]-^o a function of 

t(Lo/^o)^^'^ for experimental data. For comparison we give our results of figure 
[Hlfor Lq = 16 and 32. our numbers have been multiplied by r4He,(t>4 = 2.57. For a 
discussion see the text. 
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[C{t,Lo) — Cbuik{to)]LQ°' as a function of t{Lo/^o^2ndy^'^ , where the specific heat 
is given in units of Jmole~^i^~^ and Lq in A. To this end we have used the data 
of |l2] for thin films of "^He at vapor pressure of the thicknesses 483, 2113, 6918 
and 9869A. These data are taken from the web page [43]. For a better readabihty 
of the figure we do not give the data for 1074 and 5039A. In addition we give the 
results of [S]. In this case, we have computed Cbuikito) from the results of the fits 
given in ref. [51]. For comparison we have taken our results for Lq = 16 and 32 
from figure [TH In order to match with the experimental results, we have multiplied 
our numbers by r4,He,<i>4 = 2.57. In the high temperature phase, the experimental 
results fall nicely on top of each other. In contrast, in the low temperature phase, in 
particular for temperatures below the position of the peak, we see some scattering 
of the experimental results. In the main our results are compatible with those of 
the experiments on ^He films. In the high temperature phase our result is slightly 
larger than the experimental one. 

Finally let us compare our results with those obtained from field theory and 
from previous Monte Carlo simulations. In figure 1 (a) of [15] a one-loop result 
for the finite size scaling function fi is given. The specific heat is given in units 
of Jmo\e~^ and the length in units of A. The function has similar qualitative 
features as our result. However the peak in the low temperature phase is much more 
shallow than in our case. The maximal value is about 7, while we get 11.9. The 
position of the peak {tLy^ ^ 9) slightly differs from ours {tLy^ ^ 10.9). In [18] the 
function /i has been calculated from Monte Carlo simulations of the standard XY 
model in three dimensions. The authors have used staggered boundary conditions 
in order to suppress the order parameter at the boundary. They have simulated 
lattices of the thicknesses Lq < 24 and Li = L2 < 100. They define the specific 
heat as the second derivative of the free energy with respect to the temperature. 
They compute it using 

C = [{E') - {E)'] . (69) 

Their final result is given in figure 4 of [18] . See also, e.g., figure 2 of [Sj. From 
the scattering of the data, we conclude that statistical errors are much larger than 
in our case. Therefore the authors were not able to conduct a detailed analysis 
of corrections to scaling as we did here. Their result looks qualitatively the same 
as ours. The value of the peak seems to be somewhat larger than in our case. 
Also the position of the peak is slightly different; from figure 4 of [18j we read off 
tL^^'^ ^ —9. The authors of [19j have simulated the standard XY model using films 
of the thicknesses Lq = 12, 14 and 16. Throughout, they have used Li = L2 = 5Lq. 
They have used free ("open" in their notation) boundary conditions in the short 
direction. Their result for /i is presented in figure 2 of [19]. From the scattering 
of the results it is clear that statistical errors are much larger than in our case. 
The results of [18] and [19] are consistent. There is also reasonable agreement with 
experimental results. Similar to [IB] the value of the maximum seems to be larger 



41 



than in our case, also the position of the peak is shghtly different {tL^ ^ —8.5). 
Also in [in] corrections to scaling are not discussed. 

7 Summary and Conclusions 

We have studied the finite size scaling behaviour of the specific heat of thin films 
in the three-dimensional XY universality class. To this end we have simulated the 
improved two-component 0'^ model on the simple cubic lattice. In order the get 
a vanishing order parameter at the boundary, which is observed in experiments 
on films of ^He near the A-transition, we have employed free boundary conditions. 
These can be interpreted as Dirichlet boundary conditions with the value = (0, 0) 
of the field at the boundary. We discuss how leading boundary corrections affect 
the finite size scaling behaviour of the specific heat of thin films. We point out that 
the analytic part of the specific heat might suffer from boundary corrections that 
are not described by Lo ^ff = Lq + L^, which characterizes the leading corrections 
to the singular part. 

First we have performed simulations to get the energy density of the three- 
dimensional system for a large number of temperatures. These simulations supple- 
ment those of [27t [28] . Using these data we computed accurate estimates of the 
specific heat in the range 0.49 < (3 < 0.58 of inverse temperatures. 

Next we have analysed in detail the finite size scaling behaviour of the specific 
heat of thin films at the A-transition. To this end we have simulated films up to 
a thickness of Lq = 64 lattice units. Our result is in nice agreement with that 
obtained for thin films of ^He at the A-transition [12] . 

Furthermore we have simulated films of the thicknesses Lq = 8, 16 and 32 for a 
large range of inverse temperatures f3 in the neighbourhood of the A-transition. We 
have taken great care to obtain reliable estimates for the two-dimensional thermo- 
dynamic limit of the thin films. Using our data we have computed the finite size 
scaling functions /i and /2 defined in the introduction. It turns out that corrections 
to scaling which are caused by the free boundary conditions, have to be taken into 
account, to get a good collapse of the data obtained for different thicknesses of the 
film. These corrections can be described, to leading approximation, by an effective 
thickness -Z^o.e// = Lo + Lg of the film. In [lU] we have obtained = 1.02(7) from a 
finite size scaling study at the critical point of the three-dimensional system. How- 
ever, a priori, this i^o,e// only applies to the singular part of the specific heat. Our 
analysis of the data shows that in fact the analytic part of the specific heat requires 
an additional correction which is oc Lq^. 

The comparison of our results for the finite size scaling functions /i and /2 
and those obtained from experiments on films of ^He at the A-transition in general 
show nice agreement. We think that, in order to explain the minor discrepancies, 
a detailed knowledge of the experimental work is required. Therefore we abstain 
from any speculation on the sources of these discrepancies. 



42 



We have also compared with results obtained from field theory and previous 
Monte Carlo simulations. Field theoretic calculations are of low order; 0(e) in 
the case of the e-expansion and one or two-loop in the case of the perturbative 
expansion in three dimensions fixed. Previous Monte Carlo simulations are effected 
by relatively large statistical errors. Corrections to finite size scaling were not 
discussed in these works. 
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